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ABSTRACT  OF  THESIS 


VELOCITY  DISTRIBUTIONS  IN  THE  SEPARATED  FLOW 
BEHIND  A  WEDGE  SHAPED  MODEL  HILL 

A  portion  of  the  velocity  distributions  in  an  incompressible 
turbulent  separated  flow  behind  a  two-dimensional  model  hill  is 
investigated  experimentally.  A  2"  x  2"  wedge  is  used  as  a  model 
hill.  The  mixing  region  investigated  lies  between  the  crest  of  the 
hill  where  the  flow  separated,  and  the  reattachment  point  where  the 
separated  flow  reattaches  itself  to  the  floor. 

Gortler's  half-jet  mixing  theory  is  used  to  analyze  the 
separated  flow  which  is  curved  by  the  action  of  the  pressure  gradient 
across  the  separated  flow.  The  width  of  the  mixing  region  was  found 
to  spread  linearly  with  the  distance  as  expressed  by  the  half-jet 
theory.  If  a  suitable  similarity  parameter  is  chosen,  the  theoretical 
velocity  distribution  of  half-jet  agrees  satisfactorily  with  the  experi¬ 
ments,  except  in  the  region  near  the  floor  and  directly  behind  the 
hill.  However,  it  was  found  that  the  similarity  parameter  is  not  a 
constant,  but  is  proportional  to  the  one-half  power  of  the  distance 


downstream  from  the  hill.  Corrections  were  applied  to  account  for 


the  effects  of  the  initial  boundary  layer  thickness  and  for  the 
curvature.  The  corrections  are  based  on  Sawyer's  first  order 
theory  and  Kirk's  approximation,  but  the  results  indicate  that  these 
effects  should  be  investigated  further. 
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Chapter  I 


INTRODUCTION 

When  air  flows  along  a  flat  plate,  it  develops  a  boundary 
layer  with  a  thickness  that  increases  with  distance.  If  the  air  flow 
encounters  an  abrupt  obstacle,  such  as  a  sharp  crested  hill,  the  flow 
field  is  distorted.  The  flow  separates  at  the  crest  of  hill  and  standing 
eddies  appear  in  front  and  behind  the  hill.  The  separation  streamline 
reattaches  to  the  floor  at  the  reattachment  point  downstream  from  the 
hill.  As  shown  in  Fig.  1  ,  the  flow  field  between  the  crest  of  the  hill 
and  the  reattachment  point  can  be  divided  roughly  into  three  l  egions 
overlaying  each  other:  the  separation  bubble,  the  separation  region 
in  the  neighborhood  of  the  separation  streamline,  and  the  free  stream 
region.  Within  the  separation  bubble,  the  velocity  is  very  small  and 
may  be  considered  as  being  almost  zero.  In  the  separation  region  the 
velocity  varies  rapidly  from  zero  to  the  free  stream  velocity. 

It  is  the  objective  of  this  thesis  to  investigate  this  separation 
region  and  to  determine  if  it  is  possible  to  use  an  analysis  based  on 
tin*  half-jet  model  of  free  turbulent  shear  flow  for  the  ve!o<..it> 
distribution  in  the  separation  region.  1’he  air  flow  is  two  dimensional, 
turbulent  and  incompressible.  The  full  is  simulated  by  a  c  onstant 


2 


cross  section  wedge  2  x  2  in.  placed  perpendicular  to  the  flow 
direction  on  a  wind  tunnel  floor. 

The  separated  flow  can  be  looked  upon  as  a  free  turbulent 
shear  flow  similar  to  wakes,  jets  or  half-jets,  because  both  the 
separated  flow  and  tne  free  turbulent  shear  flow  are  similarly 
characterized  by  two  regions  with  different  velocities  and  a  mixing 
layer  within  which  the  velocity  difference  is  smoothed  out.  An 
examination  of  the  characteristic  of  each  type  oi  free  turbulent  shear 
ilow  indicates  that  the  half-jet  is  the  most  logical  model.  The  ob¬ 
jection  to  the  jet  model  can  be  seen  from  the  decrease  of  the  maximum 
jet  velocity  with  distance,  while  the  maximum  velocity  in  the  sepa¬ 
rated  flow  is  nearly  constant.  The  wake  model  is  not  desirable  either 
because  the  velocity  difference  between  the  w^.ke  and  the  free  stream 
diminishes  with  distance  which  is  not  true  for  the  separated  flow. 
Besides,  the  large  downstream  distance  required  by  available  wake 
analysis  can  not  be  obtained  in  separated  flow.  However,  the  above 
objections  vanish  when  the  half-jet  analysis  is  used.  Consequently, 
the  turbulent  half-jet  mixing  flow  io  employed  in  this  thesis  as  the 
model  of  the  separated  flow  behind  a  two-dimensional  hill  and 
Gortler's  half-jet  theory  is  applied  to  the  analysis. 

The  turbulent  half-jet  mixing  flow  has  been  studied  by  many 
investigators.  A  review  of  work  on  the  turbulent  half-jet  problem  is 
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given  in  Chapter  II.  However,  it  is  only  in  recent  years  that  the  half¬ 
jet  mixing  layer  was  employed  as  the  model  of  separated  flow  by 
Chapman  (2)  and  Korst  (7)  .  Since  then  a  considerable  amount  of  work 
concerning  the  separated  flow  in  supersonic  flow  have  been  based  on 
this  model.  Such  investigations  have  dealt  principally  with  the  base 
pressure  behind  blunt  bodies  and  the  pressure  rise  at  the  reattachment 
point  in  supersonic  flow.  Little  has  been  done  to  study  the  half-jet 
model  itself  in  the  analysis  of  separated  flow. 

With  this  in  view,  this  study  investigates  the  nature  of  the 
separated  flow  compared  to  the  half- jet  and  the  half- jet  model  is 
verified  as  usuable  for  analysis  of  the  separated  flow  region.  In 
addition,  the  flow  separating  at  the  crest  of  a  model  hill  exhibits 
certain  features  different  from  conditions  found  in  classical  half-jet. 
The  presence  of  the  lower  boundary  gives  rise  to  a  boundary  layer 
into  which  the  obstacle  is  immersed.  The  separation  phenomenon 
results  in  a  pressure  gradient  both  in  the  axial  direction  and  perpen¬ 
dicular  to  the  half- jet  flow.  An  attempt  will  be  made  to  account  for 
these  effects  by  pointing  out  possible  corrections.  A  reasonable 
consistent  model  is  developed  whose  numerical  coefficients  are 
obtained  from  the  experimental  results. 


4 


Chapter  II 

REVIEW  OF  LITERATURE 

The  purpose  of  this  chapter  is  to  present  available  information 
concerning  the  half-jet  mixing  of  turbulent,  incompressible  flow  and 
the  contributions  of  various  investigators  to  half-jet  analysis,  To 
grasp  the  full  half-jet  mixing  concept,  which  is  one  of  the  free  tur¬ 
bulent  shear  flow,  it  is  felt  that  free  turbulent  shear  flow  should  be 
discussed  first. 


A.  Free  Turbulent  Shear  Flow 

Free  turbulent  shear  Rows  occur  when  there  is  no  direct 
influence  of  a  solid  boundary.  The  effect  of  viscosity  is  very  small 
and  the  flow  spreads  laterally  into  the  ambient  fluid.  It  has  been  found 
that  free  turbulent  shear  flow  has  a  nature  similar  to  a  boundary 
layer.  This  means  that: 

1.  The  mean  flow  velocity  transverse  to  the  main  flow  is  very 
small  when  compared  with  the  main  flow  velocity. 

2.  Changes  of  quantities  in  the  direction  of  the  main  flow  are 
correspondingly  slow  with  respect  to  those  in  the  trans¬ 


verse  direction. 
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3.  The  mean  pressure  variation  across  the  flow  region  in  the 
transverse  direction  is  small. 

Consequently,  the  equation  of  motion  for  a  boundary  layer,  can  be 
used  to  analyze  a  two-dimensional  free  turbulent  flow  of  an  incom¬ 
pressible  fluid.  For  the  steady  case,  the  equation  is 


-  3u  -  8u  1  9p  _1_  9r 
dx  dy  p9x  p  3y 


(2-1) 


B. _ W.  Tollmien's  Analysis  and  Kuethe's  Extention 


Many  investigations  have  been  carried  out  in  an  attempt  to  gain 
a  knowledge  concerning  the  behavior  of  turbulent  mixing  in  an  incom¬ 
pressible  fluid.  In  1926,  W.  Tollmien  (26)  was  the  first  to  succeed  in 
the  analyzing  of  a  turbulent  mixing  flow.  He  worked  out  three 
different  types  of  mixing:  half-jet,  plane  jet  and  axially  symmetric 
jet.  The  half-jet  is  a  parallel  stream  mixing  with  the  adjacent  fluid 
at  rest.  A  plane  jet  is  formed  by  a  two-dimensional  jet  issuing  from 
a  narrow  opening  into  a  fluid  at  rest.  In  his  analysis  of  the  turbulent 
half-jet  with  zero  pressure  gradient,  W.  Tollmien  employed  the 
mixing  length  theory  and  showed  that  the  width  of  the  turbulent  jet 
increased  linearly  with  x  .  He  then  used  the  equation: 


Du 
0  x 


,  .2  Oil  ()  -  U 

Jv  oT7 


u 


♦  V 


d  U 

IK 


(2-2) 
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in  which  the  mixing  length  l  is  assumed  to  be  proportional  to  the 
distance  l  -  cx  .  c  is  a  constant  to  be  determined  from  experi¬ 
ments.  By  introducing  a  stream  function  ip  = J  u  dy  and  new  vari¬ 
ables  successively 


n *=n*  -  n*  .  n*  *  . -'-p  .  n  =  ^ 

(2c2)3 


he  obtained  the  solution  u  =  u^  F'  { r ;)  where  the  prime  denotes 
differentiation  with  respect  to  n  and 


F  .  ♦  d  e5*  /2  cosff  n*l  ♦  dje’’*  /2  siM  n* 


With  application  of  boundary  conditions 

at  the  inner  edge  r?  =  rj  u  =  Uj  F'  =  1 


(2-3) 


"  0 

9y 


f"  =  o 


v=  o 


f  =  n, 


at  the  outer  edge  n  r  h- 


u  =  0 


F'  *  0 


coefficients  in  the  solution  were 


F"  =  0 


d  *  -0.0062  d^=  0.987  d}  =  0.  577 


n*  :  0.  981 


n*2  *  -2.  04 


(2-4) 
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c  was  found  from  experiment  to  be  0.  0174  . 

Tollmien's  solution  has  been  verified  by  Liepmann  and 
Laufer  (9)  ,  who  measured  the  mean  velocity  distribution  in  a 
half-jet.  Good  agreement  was  found  between  theory  and  measure¬ 
ments.  For  plane  jet  analysis  W.  Tollmien  expressed  his  solution  as: 


(2-5) 


f  (n* )  =  e*  z  (n)  =  t'(n) 


"VT'f  jl08<Vz"+  ')  -  logV(Z  1) 

I 

This  theoretical  solution  gave  a  zero  radius  of  curvature  at  the  axis 
of  the  jet,  n*  -  0  ,  which  is  a  defect  in  Tollmien's  solution.  Later 
analysis  such  as  Gortler's  results  improved  this  point.  Forthmarm  (4) 
made  an  experimental  measurement  of  a  plane  jet  which  showed  good 
agreement  between  measured  velocity  profiles  and  Tollmien's  solu¬ 
tion,  except  near  the  axis  of  the  jet.  There  are  several  noteworthy 
assumptions  in  Tollmien's  analysis: 

1.  The  ambient  fluid  is  at  rest. 

2.  Pr&ndtle's  mixing  length  theory  io  employed. 


■■^3  tan 


-1 
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T  -  -p  u'v' 


9u  |  8u 

9y  1 9y 


the  mixing  length  i  is  assumed  constant  across  the  width 
of  mixing  region  and  proportional  to  the  distance  /  =  cx  . 

3.  The  pressure  gradient  in  the  main  stream  direction  is 

assumed  to  be  zero  =  0  . 

Ox 

4.  A  uniform  initial  velocity  distribution  before  mixing  is 
assumed,  so  that  the  initial  boundary  layer  thickness  is 
not  considered. 

5.  The  width  of  the  mixing  region  increases  linearly  with 
distance. 


6.  Similarity  of  velocity  is  assumed. 

If  any  one  of  the  above  assumptions  is  changed,  a  new  solution  will  be 
required.  Hence,  the  introduction  of  new  boundary  conditions;  such  as 
the  surrounding  fluid  being  in  motion,  or  the  application  of  a  new  tur¬ 
bulent  shearing  stress  hypothesis,  or  consideration  of  an  initial 
boundary  layer,  form  the  core  oi  subsequent  development  in  the  tur¬ 


bulent  mixing  flow. 

In  1935,  Kuethe  (8)  extended  Tollmien’s  solution  to  a  general 
case  of  the  mixing  of  two  parallel  streams  of  different  velocities.  The 
differential  equation  and  its  solutions  are  the  same  as  Tollmien’s 
half-jet  analysts.  Two  different  boundary  conditions  were  introduced 
(1)  u  u  at  n  n  ,  corresponding  to  an  ambient  fluid  with  a 

£  4 
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velocity  and  (2)  (suggested  by  Von  Karman)  =  -u->v2  • 

This  corresponds  to  the  assumption  that  no  external  forces  perpen¬ 
dicular  to  the  main  flow  are  acting  on  the  total  fluid  system.  Kuet.he 
found  that  the  effect  of  the  boundary  condition  change  from  u^  -  0 
to  u2  $  0  is  very  small.  Neither  the  mixing  region  width  nor  the 
velocity  profiles  were  affected.  The  effect  can  be  seen  only  from  a 
slight  displacement  in  the  position  of  the  edge.  In  1944,  Squire  and 
Trouncer  (25)  worked  out  the  case  of  the  axially  symmetrical  jet 
issuing  into  a  uniform  stream. 

C. _ Reichardt's  Inductive  Theory  and  Gortler's  Analysis 

One  common  feature  seen  in  all  turbulent  studies  is  that  a 
suitable  assumption  must  be  made  for  the  unknown  relation  between 
the  turbulent  shear  stress  and  the  mean  motion.  The  results  based  on 
these  assumptions  are  then  checked  by  a  final  comparison  with 
experimental  results.  It  turned  out  that  the  results  obtained  from 
each  turbulence  theory  did  not  differ  to  any  great  extent.  This 
indicates  that  the  mean  velocity  distribution  is  insensitive  to  the 
iorrn  assumed  for  the  turbulent  shear  stress,  Heichardt  (18)  in  1941, 
therefore,  suggested  an  inductive  theory.  Instead  of  beginning  with  a 
hypothesis  on  the  mixing  length  or  cnear  stress,  he  critically 
examined  the  voluminous  experimental  data  on  free  turbulent  shear 
flows,  and  discovered  that  the  velocity  profiles  under  consideration 
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could  be  approximated  very  successfully  by  the  error  function  or 
Gauss's  function.  He  also  found  empirically  that  if  there  was  a 


.  g^  & 

relation  such  as  uv  -  -A  (x)  ,  then  an  equation  describing 


the  velocity  distribution  in  the  free  turbulent  flow  can  be  obtained 


(2-6) 


Here  A  (x)  has  the  dimension  of  a  length  and  must  be  determined 

experimentally.  This  equation  for  free  turbulence  is  identical  to  the 

well  known  one-dimensional  heat  conduction  equation,  which  yields 

the  error  or  Gauss  function  as  a  solution.  Thus,  an  analogy  was 

found  between  turbulent  transfer  and  heat  conduction.  The  empirical 
,r  « 

relation  uv  =  -A  — is  then  an  expression  of  the  law  of  momentum 
transfer  which  states  that  the  flux  of  the  x -component  of  momentum, 
which  is  transferred  in  a  transverse  direction,  is  proportional  to  the 
transverse  gradient  of  momentum.  According  to  Fourier's  law  of 
heat  conduction,  the  heat  flux  is  proportional  to  the  temperature 
gradient.  ,\  (x)  is  then  analogous  to  the  heat  conductivity  and  might 
be  termed  as  a  momentum  transfer  length  As  a  consequence,  the 
differential  equation  of  heat  conduction  can  be  applied  to  turbulent 
transfer  of  momentum.  This  hypothesis  essentially  means  that  the 
exchange  coefficient  r  is  constant  across  each  section  of  the  tur¬ 
bulent  mixing  region  (14).  The  turbulent  exchange  coefficient  c  is 


ii 


defined  as  e  =  -^v'  .  Prandtl  suggested  that  when  using  this  hypothe¬ 
sis  e  =  kb  (  U  -  U  .  )  where  b  is  the  width  of  the  mixing  zone, 

and  k  is  a  dimensionless  proportional  factor.  The  new  hypothesis 
of  constant  exchange  coefficient  was  employed  by  H.  Gortler  (23)  in 
1942  to  re-examine  Tollmien's  problem.  A  much  simpler  mathe¬ 
matical  solution  to  the  half-jet  case  was  obtained  with  a  first  order 
approximation,  which  showed  good  agreement  with  the  experiments 
of  H.  Reichardt.  For  a  uniform  stream  mixing  with  a  fluid  at  rest, 
his  analysis  provided  a  velocity  distribution  expressed  by  the  error 
function 

jf*±(l+erfS)  (2-7) 

u  1  *■ 

where 

S  =  <r  ^  and  <r  is  a  constant  of  proportionality. 

Because  of  its  simplicity,  this  equation  has  been  used  widely  in  the 
analysis  of  separation  flow  (11)  and  in  the  analysis  of  base  pressure 
problems  (2)  .  The  zero  radius  of  curvature  of  velocity  profile  as 
expressed  in  Tollmien's  analysis  of  plane  jet  is  improved  when 
Eq.  (2-7)  is  employed.  However,  >.he  theoretical  mean  velocity  near 
the  edge  of  the  mixing  layer  is  higher  than  the  experimental  value. 
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D.  Dimensional  Argument 

Most  analyses  of  free  turbulent  flows  are  based  on  some 
hypothesis  for  the  mechanism  of  turbulence,  such  as  constant  mixing 
length,  or  constant  exchange  coefficient.  However,  these  assump¬ 
tions  are  not  always  in  agreement  with  measurements.  In  1947, 
Liepmann  and  Laufer  (9)  made  an  experimental  investigation  of  free 
turbulent  mixing  in  a  half- jet.  They  measured  turbulent  shearing 
stress  and  the  mean  velocity  proliles  across  the  width  of  a  mixing 
region.  Then,  the  mixing  length  and  the  exchange  coefficients  were 
computed  across  the  width.  Liepmann  and  Laufer  found  that  both  the 
mixing  length  and  the  exchange  coefficient  were  not  constant,  but 
varied  across  the  mixing  region.  Consequently,  they  concluded  that 
Tollmien  and  Gortler's  theory  were  based  on  invalid  assumptions. 

For  many  years,  there  has  been  a  growing  realization  that 
dimensional  analysis  can  furnish  practically  all  the  results  of  any 
mixing  theories.  C.  B.  Millikan  (12)  and  Squire  (24)  all  have 
expressed  this  point  of  view. 

In  1945,  Squire  has  shown  that  all  the  results  about  the  mean 
quantities  in  free  turbulent  shear  flow  can  be  deduced  by  means  of  a 
dimensional  analysis,  without  introducing  any  hypothesis  for  the  tur¬ 
bulent  mechanism.  Two  assumptions  were  made  in  his  analysis. 
First,  the  distribution  of  velocity  and  shearing  stress  are  similar  at 
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all  sections.  Second,  the  shearing  stress  is  proportional  to  the  square 
of  the  maximum  difference  in  velocity  across  a  section  and  is  inde¬ 
pendent  of  viscosity. 

Liepmann  and  Laufer  (9)  ,  in  their  investigation  of  a  half-jet, 
tried  to  describe  the  behavior  of  turbulent  mixing  by  means  of  the 
energy  and  momentum  integral  for  the  mean  and  fluctuating  motion. 
They  reached  the  same  conclusion  as  Squire,  and  once  more  demon¬ 
strated  that  a  complete  understanding  of  the  physical  aspects  of  tur¬ 
bulent  motion  is  not  necessary  for  the  analysis  of  the  free  turbulent 
shear  flow.  Dimensional  reasoning  allows  investigators  to  obtain 
the  overall  characteristic  of  the  turbulent  mixing  process. 

E. _ Torda's  Analysis  and  Kirk's  Approximation 

Previous  investigations  always  assumed  a  linear  variation  of 
+he  mixing  region  width  with  x  and  a  uniform  velocity  distribution 
in  the  stream  before  mixing  without  considering  the  influence  of  the 
upstream  boundary  layer.  Besides,  the  approaches  always  used  the 
transformation  from  the  x,  y-plane  into  an  n_plane,  where  17  is 

y 

proportional  to  ~  .  The  boundaries  of  the  mixing  region  correspond¬ 
ing  to  the  constant  values  of  r?  turned  out  to  be  straight  lines  owing 
to  the  transformation  used.  In  order  to  consider  the  effect  of  upstream 
boundary  layer,  Torda  (27)  in  1953  used  an  alternate  approach  in  his 
analysis  of  symmetric  turbulent  mixing  of  two  parallel  streams.  Both 
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streams,  which  were  separated  by  an  infinite  thin  plate  before 
mixing,  had  the  same  velocity  distribution  with  the  boundary  layer 
existing.  Torda  used  a  quartic  polynomial,  which  satisfied  all  the 
boundary  conditions,  as  an  approximation  of  the  velocity  profile 
in  the  mixing  zone. 

a  (r?2  *  l)2  -  (n2  -  2)  n2  *  (2-8) 


where 


is  a  non-dimensional  velocity  along 


the  x-axis.  Boundary  conditions  are 


y  ±6  u  =  ut 


y  =  0 


0 


Then  the  integral  equation  for  the  momentum  and  energy  equation 
was  applied  .  The  constant  exchange  coefficient  was  also  used  for 
the  turbulent  shear  stress.  With  the  application  of  boundary  con¬ 
ditions,  the  following  equations  were  obtained. 


d 

dx 


o 


n2  dy  -  Qj 


_d 

dx 


o 


Q  dy  =  0 


(2-9) 
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LA 

2  dx 


-3  ,  U1  d 

U  dy  "  T  dx 


u  dy  + 


+  kb(x)  Uj  ( 1  -a) 


0 


(2-10) 


Substituting  Eq.  (2-8)  into  Eq.  (2-9)  and  Eq.  (2-10)  yields  the 
thickness  of  the  nixing  region  as  a  function  of  x  . 


6(x)  = 


16a2 


12-11) 


where  a  is  the  function  of  x  ,  which  can  be  evaluated  graphically 
from  the  following  integration: 

a 

f  640a4  -  880a3  -  119a2+  478a  -  39  ,  k  , 

J  572(a-t)>  (16a2  -  11a  -  5)2  Qa  =  (x  +  c2> 

(2-12) 

in  which  c  and  are  constants  to  be  determined  from  experi¬ 
ments.  Thus,  from  Eq.  (2-11)  ,  the  boundaries  of  the  mixing  zone 
are  no  longer  straight  lines  but  are  curved  and  the  curvature  starts 
in  the  region  immediately  behind  the  plate.  However,  far  downstream 
the  velocity  profiles  in  the  mixing  zone  tend  to  become  similar  and 
the  boundaries  of  the  turbulent  mixing  zone  tend  to  be  straight  lint's 


as  is  expected. 
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In  1959,  Kirk  (6)  suggested  an  approximate  method  to  include 
the  effect  of  the  initial  boundary  layer  thickness  in  his  study  of  base 
pressure  in  supersonic  flow.  He  displaced  the  origin  of  the  mixing 
region  upstream  to  a  virtual  origin  and  assumed  that  the  turbulent 
mixing  started  to  develop  from  the  virtual  origin  instead  of  the  true 
origin.  He  was  then  able  to  obtain  a  similar  profile  in  the  mixing 
region. 


F. _ Sabin's  Analysis 

The  analyses  of  turbulent  mixing  discussed  above  were  all 
based  on  the  assumption  that  no  pressure  gradient  exists  in  the  main 
stream.  In  practice  this  is  not  exactly  true.  In  1963,  C.  M.  Sabin  (20) 
took  a  pressure  gradient  in  the  x-direction  into  consideration,  while 
he  studied  the  free  turbulent  mixing  between  two  parallel  streams 
with  different  velocities.  Sabin  did  not  assume  a  pressure  gradient 
across  the  mixing  region  in  the  transverse  direction.  Realizing  that 
different  hypotheses  of  turbulent  shearing  stress  always  lead  to  a 
different  form  of  solution,  Sabin  intended  to  derive  a  general 
solution  for  any  turbulent  shear  stress  hypothesis  as  well  as  for 
arbitrary  pressure  gradient  in  the  main  stream  direction.  For  this 
purpose  he  used  a  series  of  coordinate  transformations  to  absorb 
the  effect  of  the  pressure  gradient  and  the  apparent  kinematic 
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viscosity  e  .  He  began  with  the  equation  of  motion 


-  3u  .  -  8u  1  8p  ,  92u 

8x  9y  p  8x  ty7 


then  introduced  the  following  transformations: 


P  =  p(x)  "  PQ  +  7  P  u2 


P  = 


-IP  *-  P  ) 
2  v  1  z' 


x 

X=  L 


€  = 


UL 


y  = 


jL. 

UL 


U  = 


?<p.+  p2> 


2 


(2-13) 


x  =  j  e(t)  dt 
o 


A  _ 

0  =  <// 


Where  y  is  the  stream  function,  U  is  reference  velocity, 
L  is  reference  length,  Pq  is  datum  pressure  and  P  is  total 
pressure,  subscript  1  refers  to  the  upper  free  stream  and  2 
corresponds  to  the  lower  free  stream.  The  equation  of  motion  was 
finally  reduced  to  a  simplified  form 


3P 

)x 


=  (-) 
'  U 


3*P 

a, 

K* 


(2-14) 
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In  which  both  the  pressure  gradient  and  apparent  kinematic  viscosity 
disappeared.  A  solution  which  is  self- similar  was  obtained  for  the 
first  order  approximation 


2 


<s?  +  *!) 


1  + 


erf  n 


(2-15) 


where 


0 _ 

i  /£ 

2(*+  *Q)  ' 


*  =  J  J  K  +  u2)2  ^  dt 

o 

<i>  is  an  integral  constant. 

This  solution  is  independent  of  the  particular  choice  of  either  the 
pressure  gradient  or  apparent  kinematic  viscosity.  The  dimension¬ 
less  velocity  profile  in  the  mixing  region  is  found  to  be  a  function 
of  the  velocity  ratio  of  two  parallel  streams  only.  Velocity  profile 
obtained  by  this  theory  was  compared  with  those  obtained  by  the 
Gortler's  theory  for  the  case  of  a  zero  pressure  gradient.  Satis¬ 
factory  agreement  is  obtained  when  the  velocity  ratio  of  the  slow 
stream  to  the  fast  stream  is  greater  than  0.  5  .  However,  if  the 
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velocity  ratio  is  small,  lower  than  0.  4  ,  the  agreement  can  not  be 

expected,  because  the  simplifications  that  were  introduced  in  the 

process  of  comparison  are  no  longer  permissible.  An  interesting 

note  in  Sabin's  analysis  is  that  the  apparent  kinematic  viscosity  is 

insensitive  to  the  pressure  gradient  in  main  stream  direction.  Sabin 

concluded  that  Prandtl's  hypothesis  e  =  kb  (u  -  u  .  )  is 

max  min 

adequate  also  for  a  free  jet  flow  with  an  arbitrary  longitudinal 
pressure  gradient. 


rer's  First  Order  Theon 


In  1960,  Sawyer  (?.l)  conducted  an  investigation  of  a  curved 
two-dimensional  turbulent  jet  discharging  parallel  to  a  flat  plate 
at  some  distance  from  the  plate.  He  found  that  the  measured  velocity 
profiles  of  the  jet  as  it  curved  toward  the  plate,  were  nearly  sym¬ 
metrical;  although  ’here  was  a  pressure  diffei  cnee  across  the  jet. 
The  total  rate  of  jet  spread  was  almost  identical  to  that  of  a  plane 
jet.  It  seemed  ‘hat  the  curvature  had  only  a  small  effect  on  the 
velocity  distribution  in  the  jet.  In  order  to  get  a  deeper  insight  into 
this  effect.  Sawyer  {Zl)  extended  the  mixing  length  concept  to 
include  the  effect  of  the  curvature,  and  derived  a  first  order  theory 
for  turbo'  .  n‘  shearing  stress 
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r  =  pev 


i  da 

ay 


u 

mR 


(2-16) 


where  e>:  is  the  apparent  kinematic  viscosity  for  zero  curvature 
and  R  is  the  radius  of  curvature. 

The  first  order  mixing  length  theory  indicates  that  the 
curvature  has  a  considerable  effect  on  the  rates  of  entrainment,  but 
this  need  not  necessarily  be  accompanied  by  a  marked  deviation  in 
the  jet  velocity  profile  from  that  of  a  plane  jet,  since  there  exists  a 
corresponding  flow  across  the  locus  of  the  maximum  profile  velocity, 
which  balanced  the  effect  of  different  entrainment  rates. 


H,  Chapman  and  Korst's  Analysis 

Another  theory  similar  to  Torda's,  which  deals  with  the 
effects  of  the  initial  boundary  layer  on  the  mixing  within  a  half-jet, 
was  derived  by  A,  J  Chapman  and  H.  H.  Korst  in  1954  (t)  .  Their 
basic  approach  followed  the  Pai  analysis  of  jet  mixing  of  a  com¬ 
pressible  fluid.  The  one  exception  to  this  approach  was  that  they 
included  the  initial  boundary  layer.  Pai  (15)  in  1949,  proposed  a 
small  perturbation  method  to  Pnearizo  the  equation  of  motion  within 
the  mixing  region.  Thus,  the  equation  of  motion  was  simplified  to  a 


form  identical  to  that  of  the  heat  conduction  equation. 
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9u 

3?0 


o  D 


u 


(2-17) 


n  +  1 

X 

where  ?  - - n  lies  between  o  and  i  and  the  exact 

0  (n  +  1)  Ln 

value  of  n  depends  on  the  condition  of  mixing. 

Cq  is  the  empirical  constant  determined  from  experiments 
L  is  character  length 

Pai  also  proposed  a  representation  of  apparent  kinematic  viscosity  in 

n 

the  form  e  =  €0  (f  |  •  This  simplified  equation  of  motion  was 

applied  by  Chapman  and  Korst  to  the  prob'em  of  the  half- jet,  with  the 

for  -oo  <  y  <  0 

for  0  <  y  <  6 
J  o 

for  6  <  y  <  cd  (2-18) 

o 

for  y - ■>  -cd 

for  y - >  +  cd 


following  boundary  conditions. 


x  =  0 


-----  0 


=  f 


2  I 

6o  I 


X  -  X 


->  0 


u, 


u 

u, 
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where  6^  is  the  initial  boundary  layer  thickness.  They  obtained  the 
following  solution: 


u 


u 


1 


dt 


(2-19) 


in  which  the  dimensionless  variables  r]  and  17  are  written  as 

0 


>7  =  y 


\n5o  / 


n 


o 


4e  ? 


o  o 


1 

2 


If  the  thickness  of  the  initial  boundary  layer  is  very  thin  when 
compared  with  the  length  of  the  wake,  then  rjQ  -  0  .  The  solution 

reduced  to  ~  (1  +  erf  n)  (2-20) 

Cll 

This  solution  is  identical  to  thai  obtained  by  Gortler  using  the 
alternative  approach,  for  a  half-jet  with  no  initial  boundary  layer. 
Equation  (2-20)  was  adopted  in  various  studies  related  to  the  free 
turbulent  shear  flow.  Korst  (  ?)  ,  in  his  investigation  of  base  pressure 
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in  transonic  and  supersonic  flow,  applied  this  equation  to  describe 
the  motion  of  a  half-jet  mixing  and  succeeded  in  establishing  a 
criterion  for  the  reattachment  pressure  rise  in  terms  of  the  total 
pressure  developed  along  the  separation  streamline.  Mueller  (13) 
and  McDonald  (tl)  ,  also  used  this  distribution  of  velocity  in  their 
investigations  of  turbulent  separation  flow,  and  achieved  satisfactory 
results.  In  view  of  the  successful  applications  made  by  previous 
investigators,  the  same  equation  is  employed  throughout  this 
investigation  of  velocity  profiles  along  the  separation  region. 
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Chapter  III 

THEORETICAL  BACKGROUND 


This  thesis  employes  a  turbulent  half- jet  mixing  flow  as  a 
model  of  the  flow  along  the  separation  region  behind  a  two- 
dimensional  hill.  It  is  the  purpose  of  this  chapter  to  present  a 
theoretical  review  of  half- jet  analysis  following  the  work  of  Gortler. 

Suppose  two  parallel  streams  of  different  velocities,  u^ 
and  u^  ,  are  separated  by  a  thin  plate  along  the  negative  x-axis  and 
come  into  contact  at  x  =  0  (see  Fig.  2)  .  Downstream  from  this 
point  a  turbulent  mixing  region  will  occur.  The  velocity  of  upper 
stream  is  u^  which  is  larger  than  the  lower  stream  velocity  u^  . 

The  equation  of  motion  for  free  turbulent  shear  flow  in 
steady  state  is 


-  3u  |  -  du  1  9p  t  1  9t 

U  3x  V  By  p  9x  p  3} 


(3-1) 


(3-2) 
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where  b  is  the  width  of  the  mixing  region  and  k  is  the  constant. 
Eq.  (3-1)  reduces  to 


-'!rkb(5. 


(3-3) 


The  equation  of  continuity  is 


9u  f  9v 
8x  8y 


0 


(?-4) 


The  required  boundary  conditions  are 


y  =  +  oo 

u  = 

(3-5) 

y  =  -  00 

u=  s 

Introduction  of  the  stream  function  0  will  eliminate  the  continuity 
equation,  and  transform  Eq.  (3-3)  into 


80  8*0  80  8  20 

8y  3x3y  ox  dy* 


kb  (ut 


^  9V 

J2'V 


putting  ?  =  <r  ^ 

and  0  =  ^  +  u2)  J  f(S) 

where  <r  is  an  arbitrary  constant 

then  u:  ^(uJ  +  u2)  f  (?) 


(3-6) 
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Equation  (3-6)  leads  to  the  following  differential  equation  for  f(?) 


f{?)  f"(?)  +  2  —  ^ -  (?)  =  0 

X  ,  " 

Ui  +  U2 


(3-7) 


let  X  = 


U1  "  U2 


u  +  u 
i  2 


and  assume  that  the  mixing  region  width  is  proportional  to  distance 
b  =  cx  (3-8) 


where  c  is  a  proportional  constant. 

Equation  (3-7)  finally  takes  the  form  of 

f'"  (?)  +  2a f  (?)  f”  (?)  =  0 

with  ”■ =  T^rrar 

The  new  boundary  concitions  correspondingly  are 

? - •>  +oo  a  f'  (?)  =  1  +  X 

?  - — •>  -  ao  a  f'  (?)  =  1  -  X 


(3-9) 


(3-10) 


The  differential  Eq.  (3-9)  is  identical  to  the  Blasius  equation  for  a 
fin*  plate  at  zero  incicence,  with  different  boundary  conditions. 
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is  given  by  the  error  function 


(3-13) 


in  which  the  boundary  conditions  have  been  utilized.  Now  the  velocity 
distribution  across  the  mixing  region  for  the  first  order  approxima¬ 
tion  can  be  written  as 


u  = 


V  u2 


u  -  u 
i  2 

1  +  ^ - -  erf? 

U1  +  U2 


(3-14) 


If  it  is  assumed  that  u  =  0  ,  which  corresponds  to  the  case  of  a 

L 

uniform  free  stream  mixing  with  a  fluid  at  rest,  then  Eq.  (3-14) 
reduces  to 


—  =  j(l  +  erf  5) 


(3-15) 


where  ?  =  <r  ^ 
x 

The  parameter  cr  is  a  constant  which  must  be  determined  experi¬ 
mentally  10  give  the  best  agreement  between  tneory  and  measure¬ 
ments.  Eq.  (3-15)  is  the  basic  equation  for  a  half-jet  mixing  region, 
and  is  used  to  analyze  the  separated  flow  in  this  investigation.  When 
deriving  Eq.  (3-15)  it  was  assumed  that  there  was  no  initial  boundary 


29 


layer  thickness,  and  no  pressure  gradient  across  the  mixing  region 
in  the  direction  normal  to  the  main  stream. 

When  considering  the  effect  of  the  initial  boundary  layer 
thickness,  the  approximation  method  suggested  by  Kirk  (6)  can  be 
adopted  by  displacing  the  origin  of  a  half-jet  from  the  crest  of  the 
model  hill  to  a  virtual  origin  with  a  distance  xq  upstream  of  the 
original  origin.  Thus,  the  half- jet  wag;  assumed  to  have  developed 
from  the  virtual  origin  instead  of  from  the  hill  crest.  Let  be 

r 

the  similarity  parameter  for  the  case  of  the  virtual  origin,  then  the 
velocity  profiles  in  the  mixing  region  are  represented  as 


—  =  ~  1  +  erf  (a-  - ^ — 

2  o  x  +  x 

U1  L  '  ° 


(3-16) 


From  Eq.  (3-15)  and  Eq.  (3-16),  the  relation  between  a  and  <r q 
is  obtained 


cr 


x 


Equation  (3-17)  indicates  that  <t  will  increase  with  distance  x  if 
a  is  constant.  A  normal  pressure  gradient  across  the  mixing 
region  always  causes  the  flow  ucid  to  curve.  The  effect  of  this 
curvature  on  the  turbulent  shearing  stress  has  been  discussed  by 
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Sawyer  (21)  who  introduced  a  first  order  theory  which  accounts  for 
the  effect  of  this  curvature.  He  expressed  the  turbulent  shearing 
stress  as 


r  =  p  e+ 


(3-18) 


where  R  is  the  radius  of  the  curvature,  e*  is  the  apparent 
kinematic  viscosity  for  the  zero  curvature,  and  m  is  an  empirical 
constant  which  Sawyer  found  as  5  .  Equation  (3-18)  can  be  re¬ 
written  as 


t  =  p  e* 


4  Mi 

1  -  m  —  — 


u  1 

mR  - 
9u 


9u  9u 

ay =  p  6  ay 


The  effect  of  a  curvature  is  seen  in  the  apparent  kinematic  viscosity 


e=e*(l-  mQ) 


(3-19) 


u  t 

where  Q  =  —  — 

ay 


From  Eq.  (3-2)  and  Eq.  (3-8)  ,  it  yields 


e  -  kcx  (u  -  u  .  ) 
max  mm 
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where  a*  is  the  similarity  parameter  corresponding  to  conditions 

of  zero  curvature,  which  is  usually  taken  as  12  for  half- jet  mixing. 

The  effect  of  curvature  on  the  similarity  parameter  is  contained  in 

the  factor  -<  -Jr,,—  . 

•yl  -  mQ 

If  the  effects  due  to  curvature  and  initial  boundary  layer 
thickness  are  combined,  the  similarity  parameter  <r*  with  zero 
curvature  and  zero  initial  boundary  layer  thickness  will  be  modified 
by  Kq.  (3-17)  as  well  as  by  Eq.  (3-20)  .  From  Eq.  (3-20),  the 

similarity  parameter  <x  witn  curvature  but  without  initial  boundary 

o 

layer  thickness,  is  obtained. 
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<r  =  cr* 
o 


(3-21) 


which  is  again  modified  oy  Eq.  (3-1?)  to  include  the  effect  of  the 
initial  boundary  layer  thickness  expressed  in  terms  of  the  virtual 
origin  displacement  xq  . 

cr  =  —2—  -  — i- 

X  X 

,  ,  o 


■y 1  -  m  Q 


This  equation  indicates  how  the  curvature  and  the  initial  boundary 
layer  thickness  affect  the  similarity  parameter  tr  .  The  error  law. 
Equation  (3-15)  ,  for  describing  the  velocity  distributions  in  the 
separated  flow  region  shall  be  tested  and  Eq.  (3-22)  will  also  be 
checked  whether  it  suffices  to  account  for  the  effect  of  the  trans¬ 
verse  pressure  g  idient  and  of  the  initial  boundary  layer. 
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Chapter  IV 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 

Experiments  were  performed  in  a  low  speed  wind  tunnel  at 
the  Fluid  Dynamic  and  Diffusion  Laboratory  at  Colorado  State 
University.  Measured  quantities  consisted  of  mean  velocity  profiles, 
static  pressures,  turbulent  velocity  components  and  the  location  of 
reattachment.  The  following  paragraphs  will  provide  concise  infor¬ 
mation  about  the  equipment  and  procedures. 

A.  Wind  Tunnel 

Experiments  were  conducted  in  a  closed  circulating  type 

wind  tunnel  as  shown  in  Fig.  3  .  The  wind  tunnel  has  a  test  section 

with  a  usable  length  of  30  ft  ,  and  a  '  ross  sectional  area  of  6  x  6  ft2. 

Damping  screens,  located  in  the  4:  i  contraction  area,  yield  a  free 

stream  turbulence  level  —  of  less  than  l %  .  A  turbulent  boundary 

u 

CD 

layer  was  created  by  a  strip  of  l/4  in.  gravel  14  in.  long  preceded 
by  a  trip  fence  of  1/2  in.  saw  teeth.  These  obstructions  were  placed 
at  the  tunnel's  entrance.  The  hill  model,  a  2  x  2  in.  wedge  shown  in 
Fig.  4  ,  was  placed  on  the  floor  14  ft  downstream  from  the  wind 


tunnel  entrance. 


34 


Wind  speed  in  the  low  speed  wind  tunnel  was  controlled  with 
an  axial  fan,  with  16  blades,  which  is  driven  at  constant  speed  by  a 
75  hp  motor.  The  velocity  range  of  the  tunnel  is  about  5  to  75  fps 
and  is  adjusted  by  remotely  setting  the  variable  pitch  blades  of  the 
fan.  Experiments  were  carried  out  at  speeds  of  45  and  30  fps.  The 
pressure  gradient  along  the  center  line  of  the  wind  tunnel  before  the 
hill  model  was  installed  is  approximately  zero  at  30  fps.  No  pres¬ 
sure  correciions  were  made  after  the  installation  of  the  models, 
because  of  the  small  obstruction  which  the  model  offered  to  the  flow'. 
Static  pressure  along  the  floor  was  measured  by  means  of  pressure 
taps  embedded  along  the  center  line  oi  the  wind  tunnel  floor.  The 
arrangement  of  the  pressure  taps  can  be  seen  in  Fig.  3  , 

B.  Carriage 

Ver  tical  velocity  or  pressure  profiles  are  obtained  by  a  pitot 
tube  or  a  hm-wire  probe  mounted  on  a  remotely  controlled  position¬ 
ing  carriage.  Components  comprising  the  positioning  system  are  a 
28  V.  D.C.  motor,  f.  potentiometer,  a  screw  drive,  a  probe  holder, 
a  guide  bar  and  a  control  box  which  is  used  to  control  the  speed  and 
direction  of  the  motor.  The  D.  C.  motor  drives  the  screw  drive, 
which  in  turn  moves  the  orobc  up  and  down.  A  probe  is  attached  to 
the  guide  bar  by  means  of  probe  holder  which  allows  the  probe  to 
be  moved  a  distance  of  19  inches.  A  potentiometer  is  fixed  on  the 
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probe  holder  and  geared  to  the  guide  bar.  As  the  probe  holder 
moves,  a  voltage  drop  corresponding  to  the  distance  travelled  by 
the  probe  holder  is  produced  across  the  potentiometer.  This  voltage 
is  then  transmitted  to  the  x-axis  of  a  X-Y  plotter  to  give  the 
position  of  the  pitot  tube. 

C.  Pitot  Tube,  X-Y  Plotter,  and  Electronic  Manometer 

A  pitot  tube  was  used  to  measure  mean  flow  velocity.  The 
hot-wire  anenometer,  which  will  be  discussed  later,  is  used  to 
measure  the  fluctuating  components  of  velocity.  Outside  diameter 
of  the  pitot  tube  was  3/ 1 6  in.  Plastic  tubing  was  used  to  connect  the 
pitot  tube  to  a  Transonic  Equibar  Type  120  electronic  manometer, 
which  has  8  ranges  of  pressure  from  0.  01  mm  Hg  to  30  mm  Hg 
with  a  dc -output  of  30  mv  for  each  full  scale  reading.  The  electric 
manometer  can  detect  pressure  changes  of  0.  0002  mm  Hg.  Output 
of  the  electronic  manometer  was  transmitted  to  the  y-axis  of  a 
X-Y  ptotter,  Moseley  Type  135  .  A  voltage  corresponding  to  the 
position  of  pitot  tube  was  applied  to  the  x-axis  of  the  X-Y  plotter. 

In  this  manner  it  was  possible  to  plot  out  a  continuous  profile  of 
the  dynamic  pressure  or  static  pressure  directly. 
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D.  Mean  Velocity  Measurement 

Mean  velocity  was  measured  with  a  pitot  tube,  which  was 
mounted  on  the  carriage.  The  free  stream  velocity  was  set  at  45 
and  30  fps.  It  was  measured  by  placing  the  pitot  tube  2  ft  above  the 
floor,  upstream  of  the  hill  model.  Vertical  velocity  profiles  were 
measured  from  the  floor  to  a  height  of  19  inches  and  from  the 
model's  crest  to  a  distance  of  60  inches  downstream.  Measurements 
were  taken  at  various  stations  behind  the  model  along  the  tunnel's 
center  line.  Air  flow  two-dimensionality  was  verified  by  velocity 
profiles  taken  at  a  few  stations  located  on  either  side  of  the  center 
line.  These  profiles  were  within  one  percent  agreement  as  shown 
in  Fig.  5  .  Air  temperatures  and  barometric  pressure  were  also 
recorded  during  test  runs. 

The  mean  velocity  measured  with  a  pitot  tube  is  calculated 

from  the  corresponding  reading  of  the  transonic  manometer  from 

the  relation  u  =  2.  36  \  ■—  .  h  is  the  transonic  reading  (in  mm  Hg) 

V  P 

corresponding  to  the  dynamic  head,  p  is  the  density  of  air  (  in 


lb  -  sec 
fpr 


)  and  u  is  the  velocity  (  in  fps  )  . 


A  dynamic  pleasure  is  defined  aa  the  difference  of  the  total 


pressure  and  the  static  pressure  at  the  same  point.  However,  for 


the  pitot  tube  used,  the  ?ota!  pressure  and  the  static  pressure  arc 
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not  measured  at  the  same  point,  but  at  points  one-inch  apart  as 
shown  in  Fig.  6.  If  H  ,  represents  the  total  pressure  at  point  i  , 
and  if  p^  and  p^  represent  the  static  pressure  at  point  1  and 
point  2  respectively,  then  the  velocity  at  point  i  ,  in  terms  of 
dynamic  pressure  should  be: 

i p  flZ  =  Hi  -  Pi 


However,  the  dynamic  pressure  measured  with  the  pitot  tube  is 


1  —  2  tJ 
-  p  u  =  H  -  p 
2  m  1  2 


where  u  corresponding  to  the  velocity  measured  with  the  pitot 
tube.  The  theoretical  value  and  the  measured  value  can  be  related  as 


Ht  -  p2+  P2-P, 


=  —  p  u2  +  ;  ^  Ax 
2H  m  •  3x  1 


The  error  of  velocity  due  to  pitot  tube  measurement  is 


-  u  -  u  u 

Au _ m  .  j  _  in 


u 


u 


u 
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Ax 


If  the  dynamic  pressure  is  expressed  in  terms  of  the  pressure 
transducer  reading  h  ,  then  the  error  takes  the  form 


Au 

u 


(4-1) 


where  y  is  the  specific  weight  of  air,  and 
mated  by  the  pressure  difference,  i,  e. 


9p 

Bx 


I  1 


Ax 


Ax  is  approxx- 


there  is  no  pressure  gradient  in  the  x-direction,  the  pitot  tube  gives 
the  true  value  of  velocity.  Further,  if  the  difference  of  pressure 
Ap  is  small  when  compared  with  pitot  tube  reading  h  ,  the  error 
would  be  very  small  and  the  velocity  calculated  directly  from  pitot 
tube  reading  is  still  within  the  accuracy  required.  In  this  investi¬ 
gation,  the  maximum  error  of  velocity  with  the  free  stream  velocity 
as  reference  is  within  1.5%,  so  no  effort  was  made  to  correct  the 


pitot  tube  readings. 
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Taepmann  and  Laufer  (9)  have  used  both  pitot  tube  and  the 
hot-wire  to  measure  the  same  mean  velocity  distribution  separately, 
to  see  if  there  was  any  discrepancy  existing  between  the  results. 
They  found  both  methods  gave  satisfactory  agreement. 

E.  Static  Pressure  Measurement 

Only  the  static  hole  of  the  pitot  tube  was  used  when  measur¬ 
ing  static  pressure.  Measuring  instruments  (seen  in  Fig.  7)  were 
arranged  in  the  same  manner  as  for  the  measurements  of  the  veloc¬ 
ity  measuring  device,  except  that  another  pitot  tube  was  used  to 
give  the  static  reference  pressure.  Throughout  the  experiment,  the 
reference  pitot  tube  was  fastened  to  a  slender  bar  located  upstream 
from  the  hill  model.  The  ba  was  suspended  from  the  tunnel's 
ceiling  and  extended  approximately  20  in.  into  the  ambient  stream. 
Static  pressures  were  measured  around  me  hill,  along  the  center 
line  of  the  floor.  Pressures  around  the  hill  were  measured  by 
means  of  777  in.  pressure  taps  around  the  model,  (see  Fig.  4) 

Static  pressures  along  the  floor  were  measured  through  pressure 
taps  embedded  in  the  floor.  These  pressure  taps  were  connected 
successively  through  plastic  tubing  to  the  electronic  manometer, 
where  they  were  measured  against  the  reference  static  pressure. 
Vertical  static  pressure  distributions  were  measured  in  the  same 
w ay  as  the  mean  velocity  profiles,  except  that  only  the  static  hole 
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was  used.  Some  fluctuation  of  static  pressure  was  observed  through¬ 
out  the  measurements.  The  largest  amounts  of  fluctuation  were 
found  in  the  mixing  region. 


F.  Turbulent  Measurement 


Turbulent  quantities  u'z  ,  u'v'  and  v'  were 
evaluated  by  using  a  constant  temperature  hot-wire  anemometer. 
The  fluctuating  velocity  component  in  the  main  stream  direction 
u,z  was  measured  with  a  single  wire  probe.  A  cross  wire  was 
used  in  place  of  the  single  probe  when  measuring  u'v'  and  v' 2  . 
Turbulent  quantities  were  measured  with  the  probe  mounted  on  the 
positioning  carriage.  A  platinum  wire;  which  had  a  diameter  of 
0.  0001  in.  ,  a  length  of  0.  1  in.  and  a  cold  resistance  of  approxi¬ 
mately  5  ohms;  was  used  as  a  sensor.  This  wire  was  mounted  on  a 
hot-wire  probe  made  by  DISA  Company.  The  hot-wire  probe  was 
operated  by  a  constant  temperature  type  servo  amplifier,  Hubbard 
3A  ,  from  which  the  A.C.  output  was  led  into  a  true  rms-meter. 
Type  BRUEL  and  KJAER  2409  .  This  instrument  had  been  modified 
to  give  an  output  voltage  proportional  to  the  meter  reading.  Output 
voltage  from  the  rms-meter  and  the  voltage  proportional  to  the 
probe  position  were  then  applied  to  the  two  axes  of  an  X-Y  plotter, 
thus  yielding  continuous  turbulent  profiles.  Instrument  arrangement 
is  shown  in  Fig.  8  .  The  output  voltage  of  the  Hubbard  amplifier  is 
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a  linear  function  of  velocity.  From  measurements  with  pitot  tube 
and  hot-wire,  the  free  stream  velocity  is  calibrated  against  the 
output  reading  of  the  hot-wire  anemometer.  A  straight  line  cali¬ 
bration  curve  was  obtained.  Fig.  9  gives  an  example  of  this  cali¬ 
bration  curve.  With  the  slope  of  the  calibration  curve  determined, 
the  turbulent  quantities  are  obtained  by  conversion  from  the  reading 
of  the  plotted  profiles. 

Turbulent  velocity  measurement  made  with  a  hot-wire 
anemometer  are  based  on  the  principle  that  the  heat  loss  per  unit 
time  is  equal  to  the  heat  generated  per  unit  time  by  an  electric  cur¬ 
rent  passing  through  the  wire.  Heat  loss  is  generally  heat  conduc¬ 
tion,  radiation,  free  and  forced  convection.  In  nractice  the  effects 
of  radiation  and  free  convection  are  negligible.  By  King's  law,  the 
empirical  relation  between  the  heat  loss  and  the  mass  flow  can  be 
written  as 

I2R  __ 

rrR=3tD^  <4-2> 

w  a 

Where  R  is  the  wire  resistance  while  in  operation,  R  is  the 

W  SI 

wire  resistance  at  air  temperature  and  B  and  D  are  constants. 

For  constant  temperature  operation,  the  resistance  is  kept 
constant,  the  fluctuation  of  the  velocity  component  will  give  rise  to 
a  fluctuation  of  heating  current.  If  fluctuations  of  velocity  and 
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current  are  both  small  compared  with  the  mean  values,  then 


2  I  R 

w  ,T  _  _ 

R  -  R  dI '  2 

w  a 


dl  =  I' 
du  =  u’ 


The  relation  between  the  fluctuation  of  velocity  and  that  of  current 
is  then  determined  as 

I'  *  Su'  (4-3) 

where  S  is  the  sensitivity  of  the  hot-wire  anemometer.  In  the 
instrument  used,  Hubbard  Type  3A  ,  the  current  (in  milliampere) 
is  passing  through  a  resistor  of  50,  000  ohms,  across  which  the 
output  voltage  is  measured.  The  output  voltage  e'  (in  volts)  is 
then  related  to  the  fluctuating  current  through 

o'  =  50  I'  (4-4) 

and  to  the  fluctuating  velocity 
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where  e'  is  the  mean  square  reading  of  the  voltage  output,  and 

A  =  Jr  =  .  The  inverse  of  the  sensitivity  S  is  defined  as  the 

S  dl 

calibration  constant  A  ,  which  is  obtained  from  the  calibration 
curve  of  the  mean  velocity  plotted  against  the  mean  current.  A 
linearing  circuit  designed  into  the  instrument  is  responsible  for  the 
straight  line  calibration  curve  in  Fig.  9  .  A  is  then  determined  iy 
measuring  the  slope  of  the  calibration  curve.  As  a  result,  the  tur¬ 
bulent  quantities  u'z  can  be  obtained  from  Eq.  (4-5)  . 

Turbulent  quantities  u’v'  and  v'2  were  measured  with  a 
cross  wire  probe.  The  two  wires  comprising  this  probe  were  in¬ 
clined  at  angles  a  and  3  to  the  main  stream  direction  as  indicated 
in  Fig.  10  .  For  this  wire  arrangement,  the  relations  between 
fluctuations  of  currents  and  of  velocities  become 


AjI'j  =  u'  cos  a  +  v'  sin  a 
A -  u'  cos  3  -  v'  sin  3 

Square  and  then  take  the  time-average 


for  wire  1 


for  wire  2 


u' 2  cos4  or  +  2u'v'  sin  a  cos  a  *  v'*  sin4  a 


(4-6) 
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A  2 

■** 2  |  _ 

gO~  j  =  u'2  cos2$  -  Z  u'v'  sin8  cos8  +  v'2  sin20 

(4-7) 

For  the  experiments  presented  in  this  study  the  angles  were  a  =  38°, 
8  =  43°  and  the  calibration  curve  was  given  by  the  equation 
u  =  43  I  i  7.  5  ,  where  =  A^  =  43  .  With  these  values,  sub¬ 
tracting  Eq.  (4-7)  from  Eq.  (4-6)  yields 

777  (  el”2'  -  e\r)  =  u*^~(cos2a  -  cos28)  +  uV  (sin2a  +  sin2(3)  + 

OU  1  L 

+  v'2  (  sin2ar  -  sin2/3) 

which  yields  u'v'  =  0.  376  (  e’^2  -  e^2)  (4-8) 

Further,  adding  Eq.  (*>-<”)  and  Eq.  (4-7)  leads  to 

(  e'2+  e'2)  =  u'2  (  cos?ar  +  cos23)  +  u'v'  (sinZcr  -  sin23)  + 
|d0  1  Z 

+  v' 2  (sin2a  +  sin28) 

or  =  0.377  -  t.  37  u77"  (4-0) 

1  w 

From  Eq.  (4-3)  and  Eq.(4-9)  ,  v77'  and  u’v!  can  be  calculated. 
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G.  Measurement  of  Reattachment  Point 

A  special  technique  was  developed  (16)  to  measure  the 
reattachment  point.  Two  hot  wires  of  exactly  the  same  length  and 
resistance  were  placed  parallel  and  as  close  as  possible  without 
touching  in  the  separated  region.  Each  one  was  operated  by  one 
channel  of  the  two-channel  Hubbard  hot-wire  anemometer  amplifier 
The  outputs  of  the  hot-wire  amplifiers  were  fed  into  a  special  dis¬ 
criminating  circuit,  designed  by  Mr.  Calvin  Finn.  This  discrimi¬ 
nating  circuit  provides  a  train  of  100,  000  pulses  per  second  as  long 
as  the  input  to  one  preset  discriminator  channel  (channel  A)  is 
larger  than  the  input  to  the  other  channel  (channel  B)  .  If  the  input 
to  channel  B  is  larger  then  no  signal  is  generated  by  the  circuit. 

Now  if  the  two-wire  probe  is  placed  at  some  suitable  distance 
above  the  floor  with  the  >ia;ie  through  the  axes  of  the  wires  parallel 
to  the  floor,  then,  when  the  flow  velocity  is  parallel  to  the  floor,  the 
downstream  wire  is  in  the  wake  of  the  upstream  wire  and  conse¬ 
quently  does  not  get  cooled  as  much  as  the  upstream  wire.  Under 
these  conditions  the  fee  'back  amplifier  would  give  a  iarger  current 
output  for  the  wire  located  upstream  than  for  the  downstream  wire. 

If  the  signals  from  the  two  .ires  are  fed  into  the  discriminating 
circuit  it  yields  an  output  which  when  transferred  to  a  pulse  counter 
(Hewlett  Packard  Type  5 ZZ)  .  indicates  the  percent  of  tune  one  wire 
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yields  a  larger  signal  than  the  other.  The  reattachment  point  was 
defined  as  that  location  where  the  individual  direction  of  flow  is 
equal  in  both  directions  {  a  50%  reading  on  the  discriminator)  .  The 
reattachment  point  was  thus  found  at  25.  5  in.  downstream  from  the 
hill  when  the  free  stream  elocity  was  45  fps.  The  equipment  is 
shown  in  Fig.  1 1  . 


4  7 


Chapter  V 

ANALYSIS  AND  DISCUSSION 

The  flow  along  the  separation  region  behind  a  hill  is  similar 
to  a  half-jet  mixing  with  the  ambient  fluid  at  rest.  Therefore,  the 
velocity  profiles  along  the  separation  region  can  be  described  by 
Equation  (3-15)  which  was  obtained  for  a  half-jet 

—  *  {(1  +  erf  I) 

U1 

y 

where  ?  =  <r  and  <r  is  a  constant, 
x 

In  this  equation,  x  and  y  are  the  coordinates  of  an  intrinsic 
coordinate  system  in  the  flow.  The  relation  between  the  intrinsic 
coordinate  and  the  reference  coordinate  system  can  be  derived 
from  Eq.  (3-15)  ,  as  will  be  discussed  below. 

A.  Free  Stream  Velocity 

When  air  flows  over  the  hill,  a  separation  bubble  forms 
directly  behind  the  hill.  Inside  the  bubble  there  exists  a  vertex 
region  where  the  velocity  is  very  small  and  may  be  considered  as 
being  almost  equal  to  zero.  Between  the  bubble  and  the  free  stream 
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there  is  a  mixing  layer  across  which  the  velocity  increases  rapidly 
from  zero  to  the  free  stream  velocity.  Influence  of  the  hill  causes 
the  velocity  in  the  separation  region  above  the  bubble  to  rise  to  a 
peak  value  near  the  crest  of  the  hill,  and  then  decrease  gradually 
over  a  distance  to  reach  an  equilibrium  value.  The  amount  of  this 
velocity  decrease  is  about  5%  to  8%  of  the  maximum  value.  Maxi¬ 
mum  velocity  is  then  defined  as  free  stream  velocity  u^  .  Mean 
velocity  profiles  behind  the  hill  are  shown  in  Fig.  12  .  The  veloci¬ 
ties  near  the  floor  are  left  out  because  there  is  a  reversed  flow 
near  the  floor,  which  is  hard  to  measure  correctly  with  the  pitot 
tube.  The  locus  of  the  points  which  correspond  to  the  maximum 
velocity  is  then  defined  as  the  upper  boundary  of  the  half- jet  mixing 
region,  which  is  a  small  gradient  curve  seen  in  Fig.  13  .  With  the 
free  stream  velocity  determined,  the  ratio  of  local  velocity  to  the 
free  stream  velocity  is  only  a  function  of  local  position  as  described 
by  Eq.  (3-15)  . 

Figure  14  presents  the  velocity  variation  along  the  upper 
boundary  where  ^  =  1  ,  and  along  the  central  velocity  line 

—  =  j  ,  that  is,  the  jet-axis.  The  velocity  was  made  dimension- 

U1  Z 

less  by  dividing  by  the  corresponding  velocity  at  the  crest  of  the 
hill.  A  small  variation  of  the  half-jet  free  stream  velocity  was 
observed  along  the  distance  downstream.  At  first,  the  velocity 
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increases,  then  it  decreases  and  tends  to  approach  a  constant  vaiue 
downstream  near  the  end  of  the  separation  region.  The  pressure 
gradient  in  the  x-direction,  which  is  implied  in  the  variation  of  the 
free  stream  velocity,  was  shown  by  Sabin  (20)  as  making  only  a 
small  contribution  to  the  turbulent  shear  stress.  If  the  variation  of 
free  stream  velocity  is  neglected,  which  means  the  pressure  gra¬ 
dient  in  the  x-direction  is  assumed  to  be  zero,  the  analysis  is 
simplified  and  the  error  introduced  would  not  have  a  significant 
effect  on  the  end  results.  This  has  been  done  in  this  analysis. 

B.  Transformation  of  Coordinates 

An  intrinsic  coordinate  system,  x  and  y  ,  is  used  in  half- 

j 

i 

jet  analysis,  refer  to  Eq.  (3-15)  .  Measurements  use  a  fixed 
reference  coordinate  system.  This  system  consists  of  X  and  Y-axis 
which  originated  at  the  front  edge  of  the  hill  model.  The  X-axis 
appears  along  the  floor  with  a  positive  downstream  direction.  The 
Y-axis  is  normal  to  the  floor  with  a  positive  direction  upward.  The 
relation  between  the  coordinate  system  can  be  obtained  from  Eq. 
(3-15)  . 

Location  of  the  X-axis  can  be  found  from  the  experimental 
data  according  to  the  following  procedure.  If  the  velocity  distribution 
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is  given  by 

=-  =  i(l  +  erf?) 

U1 

where  $  =  o-  ^ 
x 

then  it  follows  that 

if  y  =  0  then  ?  =  0 

and  the  velocity  ratio  becomes 

u  _  _ 

*1  ~  2 

This  means  that  the  x-axis  coincides  with  central  velocity  line 

where  ~  ~  ^  •  Since  the  free  stream  velocity  has  been 

U1 

determined  as  the  maximum  velocity,  the  x-axis  can  be  found  easily 
from  the  measured  velocity  profile  by  locating  the  points  where 

—  =  .  Fig.  13  illustrates  an  x-axis  obtained  in  the  foregoing 

ui 

manner.  This  x-axis  was  observed  to  curve  up  from  the  crest  of 
the  hill  and  soon  approach  a  straight  line  parallel  to  the  floor.  The 
angle  between  the  y-axis  and  the  Y-axis  due  to  the  curvature  of  the 
x-axis  is  so  small  as  to  be  negligible.  Therefore,  the  y-axis  is 
assumed  to  coincide  approximately  with  the  Y-axis.  Let  Yc 
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represent  the  Y -coordinate  of  the  central  velocity  line,  namely,  the 
x-axis,  then  a  relation  between  y  and  Y  can  be  established  as 
y  =  Y  -  Y^  shown  in  Fig.  15  . 

The  initial  curvature  section  of  the  x-axis  is  rather  short. 
Thus,  the  distance  along  the  x-axis  is  approximated  by  the  distance 
along  the  X-axis,  i.  e. ,  by  the  distance  parallel  to  the  floor.  The 
error  introduced  by  such  an  approximation  is  small  at  large  distance. 
Finally,  the  transformation  of  the  intrinsic  coordinate  system  into 
the  fixed  reference  system  is  given  by 

y  =  Y  -  Yc  (5-1) 

x  =  X 

At  this  point,  it  is  possible  to  compare  the  measured  value  with  the 
theoretical  value  obtained  by  Eq.  (3-15)  ,  if  there  is  a  suitable  value 
for  similarity  parameter  <r  . 

C  Half-Jet  Boundary 

Theoretically,  the  velocity  profiles  given  by  Eq.  (3-15) 
result  in  a  velocity  equal  to  the  free  stream  velocity  at  an  infinite 
distance  from  the  x-axis.  Therefore,  the  boundary  line  is  usually 

defined  at  =  0.99  to  avoid  this  difficulty  In  this  investigation, 

^  4 


the  upper  boundary  is  defined  as  the  maximum  velocity  point  in  the 
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profile.  The  upper  boundary  is  presented  in  Fig.  13  ,  which  does 
not  appear  as  a  straight  line  but  as  a  curve.  This  curvature  is  large 
directly  behind  the  hill,  but  at  a  short  distance  downstream  the 
curvature  decreases  to  the  point  that  it  may  be  approximated  by  a 
straight  line.  This  curved  boundary  is  in  contradiction  to  the 
classical  half-jet  studies.  However,  in  an  investigation  of  turbulent, 
incompressible,  symmetrical  mixing  of  two  parallel  streams, 

T.  P.  Torda  also  found  a  curved  boundary  which  he  attributed  to 
the  initial  boundary  layer  thickness  upstream  of  the  trailing  edge  of 
the  flat  plate.  In  addition  to  the  effect  of  initial  boundary  layer 
thickness,  the  obstruction  of  the  hill  in  the  flow  might  contribute  to 
the  curvature  of  the  half- jet.  The  hill  extended  from  the  floor  into 
the  flow.  Therefore,  the  flow  upstream,  below  the  height  of  the  hill, 
was  forced  to  rise  as  it  approached  the  hill.  As  a  consequence,  a 
velocity  component  normal  to  the  mean  stream  direction  appeared 
near  the  hill.  The  half-jet  mixing  region  which  was  slightly  displaced 
by  this  vertical  velocity  component  will  have  curved  appearance.  This 
normal  velocity  died  out  further  downstream  so  that  the  curvature 
reduced  correspondingly. 

The  lower  boundary  for  half- jet  mixing,  which  lies  within 
the  vortex  region  in  the  separation  bubble,  was  not  well  defined.  An 
attempt  was  made  to  define  this  lower  boundary  by  a  constant  velocity 
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ratio  line  3—  =  0.  1  ,  but  the  data  are  too  irregular  to  define  a 
Ui 

regular  curve.  There  is  a  reversal  of  flow  within  the  separation 
bubble  caused  by  the  entrainment  of  air  into  the  jet  region.  The 
pitot  tube  was  not  oriented  to  measure  this  reverse  velocity.  The 
large  scatters  in  the  reading  may  be  due  to  the  inability  of  pitot  tube 
to  read  accurately  the  velocity  in  the  bubble. 

D.  Spread  of  Width 

In  view  of  the  difficulty  in  defining  the  lower  boundary  of  a 
half-jet  mixing  region  as  mentioned  before,  the  mixing  region  is 
divided  into  two  zones:  Zone  I  and  Zone  II  .  Zone  I, as  seen  in  Fig.  16, 
is  the  region  above  the  x-axis  up  to  the  upper  boundary  line  where 

^—  =  i;  and  Zone  II  is  the  region  below  the  x-axis.  Both  the  upper 
U1 

boundary  and  the  x-axis  of  the  half-jet  are  curved  rather  than  linear. 
The  distance  between  the  twolines,  which  is  the  width  of  Zone  I  ,  was 
observed  to  vary  approximately  linearly  over  a  distance  of  4  to  20 
inches  downstream  from  the  model  hill.  Further  downstream,  it 
tends  to  approach  an  asymptotic  value.  This  linear  variation  of  the 
width  with  distance  is  in  agreement  with  half- jet  theory. 

The  spread  of  the  width  of  Zone  I  is  plotted  in  Fig.  17.  The 
width  of  Zone  II  is  no*  shown  because  the  low  boundary  of  the  jet  is 
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poorly  defined.  The  spread  of  the  width  of  the  mixing  region  is  found 
to  decrease  with  an  increase  of  free  stream  velocity.  The  width 
when  free  stream  velocity  is  30  fps  is  bigger  than  when  the  velocity 
is  45  fps. 


E.  Proximity  of  Floor 

In  half-jet  analysis,  it  is  supposed  that  no  solid  boundary 
exists  in  the  neighborhood  of  the  mixing  region.  However,  in  this 
study  of  separation  flow,  there  is  a  floor  behind  the  model  hill.  The 
effect  due  to  the  proximity  of  the  floor  is  felt  in  the  region  near  the 
floor  and  will  decrease  with  the  distance  from  the  floor. 

Zone  I  is  far  enough  away  from  the  floor  that  the  effect  of 
the  floor  is  so  slight  as  to  be  unnoticeable,  and  the  separation  flow 
in  this  region  develops  as  a  half-jet.  In  Zone  II  ,  which  is  below 
the  jet-axis  and  near  the  floor,  the  influence  of  the  floor  on  the  flow 
is  stronger.  This  influence  is  found  from  measured  velocity  data 
deviation  from  analytical  results.  The  floor  also  restricts  the  width 
of  Zone  II  and  as  a  result  Zone  II  is  narrower  than  Zone  I  . 

The  most  important  floor  effect  is  the  formation  of  a  separa¬ 
tion  bubble  which  adds  complexity  to  the  problem.  Jet  flow  near  the 
outer  edge  of  the  bubble  entrains  the  fluid  from  the  bubble  by  the 
action  of  the  shearing  stress,  thus  causing  the  pressure  inside  the 
bubble  to  decrease.  This  decrease  in  pressure  forces  the  lower  part 
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of  half-jet  toward  the  floor.  The  bubble  also  causes  the  half-jet 
mixing  region  to  spread  vertically  as  distance  increases.  As  a 
result  of  this,  there  is  a  point  downstream  from  the  hill  where  the 
flow  reattaches  to  the  floor  with  an  amount  of  the  flow  reversed 
back  into  the  bubble.  Mass  conservation  requires  that  the  amount  of 
fluid  entrained  from  the  bubble  should  be  balanced  by  the  amount  of 
the  reversed  flow  at  reattachment.  The  air  in  the  separation  bubble 
is  thus  induced  into  a  circulatory  motion  to  form  an  eddy  vortex. 
Because  it  is  hard  to  use  pitot  tube  to  measure  the  reversed  flow  in 
the  bubble,  the  velocities  near  the  floor  are  left  out  in  Fig.  5  and 
Figure  12  .  However,  the  pitot  tube  is  capable  of  measuring  the 
static  pressure  in  the  reversed  flow,  and  the  pressure  inside  the 
bubble  is  found  lower  than  the  pressure  outside. 

F.  Pressure  Gradient  and  Curvature 

Pressure  gradients  are  found  both  in  the  main  stream  direc¬ 
tion  and  in  the  transverse  direction.  Fig.  18  shows  the  variation  of 
static  pressure  along  the  constant  velocity  ratio  lines  and  Fig.  19 
shows  the  variation  of  static  pressure  at  constant  heights  above  the 
floor,  downstream  of  the  hill.  The  development  of  static  pressure 
distribution  downstream  of  the  hill  is  shown  in  Fig.  20  .  Near 
reattachment  the  static  pressure  was  observed  to  increase  over 
a  large  horizontal  distance.  Below  the  jet-axis  the  pressure  is 
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nearly  uniform.  Most  of  the  variation  of  static  pressure  from  the 
pressure  in  the  separation  bubble  to  the  pressure  in  the  external 
free  stream  takes  place  in  the  region  above  the  jet  axis  to  the  outer 
jet  boundary,  that  is,  Zone  I  .  The  pressure  difference  across 
Zone  I  becomes  smaller  and  smaller  further  downstream.  The 
horizontal  pressure  gradient  has  been  found  by  Sabin  (20)  to  have 
only  a  slight  effect  on  the  turbulent  mixing  mechanism  and  conse¬ 
quently  was  neglected  in  this  investigation.  However,  transverse 
pressure  gradient  across  the  mixing  region  induces  a  curvature  in 
the  separated  flow. 

The  curvature  can  be  seen  in  the  equation  in  the  cylinder 
coordinates  (5)  .  Assume  u  is  parallel  to  the  streamlines,  so  that 
v  =  0  .  Then  the  equation  of  motion  along  a  streamline  is 

u 2  1  op  .  ,  v73"  ,  1  9uV  urr  ,x 

—  - i_  -t-  -  +  -  + - - -  (a- 2) 

r  p  9r  9r  r  r  90  r 

in  which  the  viscous  stress  components,  which  are  very  small  in 
free  turbulent  shear  flow,  are  omitted.  Eq.  (5-2)  may  be  rearranged 
as 

ui  r  u v' *  _  1  0(u' v ')  ]  1  9p  t  9 v 1  2 

r  [  u*  u2'  ft7  30  J  p  9r  3r 

Assume  the  streamlines  in  the  separated  flow  are  parallel  to  each 
other,  so  the  radius  of  streamline  K  is  constant  at  each  cross 
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section,  i.  e. ,  R  is  not  a  function  of  r  .  Integrating  the  above 
equation  over  r  , 


R 


9u' v1 
9(9 


dr  = 


J 


9p  ,  s'v*" 

4.  p  - - 

9r  p  9r 


dr 


u  v 

As  shown  in  Fie.  21  ,  the  turbulent  quantities  zri —  ,  =rr~  ,  ~ 

49  u  u  u 


u'v' 


are  very  small.  Their  integral  are  much  smaller  when  compared 
with  the  integral  of  the  ^rm  u2  ,  so  these  turbulent  quantities  are 


neglected.  Ii  yields 


(5-3) 


The  radial  coordinate  can  be  approximated  by  the  y-coordinate.  The 
integral  is  taken  from  the  jet-axis  to  the  outer  jet  boundary.  Then 
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Since  p  Av^2”  is  much  smaller  than  Ap  as  can  be  seen  in  Fig.  21, 
it  is  neglected  .  .  The  radius  of  curvature  is  thus  obtained 


(5-4) 


The  pressure  difference  Ap  decreases  with  x,  so  the  radius  of 
curvature  R  increases  with  x  as  shown  in  Fig.  22  .  This 
relation  indicates  how  the  separated  flow  is  bent  by  the  action  of 
the  transverse  pressure  gradient. 


G.  Similarity  Parameter 

The  ernperical  similarity  parameter  or  is  usually  chosen 
when  fitting  velocity  data  to  theory.  Previous  investigators  have 
found  various  values  of  <r  .  For  instance,  Liepmann  and  Laufer  (9)  , 
in  their  experiments  of  a  two-dimensional  mixing  layer,  found  that 
the  best  value  for  presenting  experimental  results  according  to 
Tolimien'j  theury  was  12  .  but  by  Gortler's  theory  <r  =  1 1  gave 
better  agreements.  Reichardt  (18)  found  that  a  =  13.5  for  the 
half-jet  case  and  cr  =  7.  67  for  the  two-dimensional  jet  when 
Gortler's  theory  was  employed.  Sawyer  (2*  cl.ose  c r  =  15  in  his 
analysis  of  a  two-dimensional  jet  issuing  parallel  \o  a  flat  plate, 
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The  non-dimensional  velocity  profiles  are  given  by  Eq.  (3-15) 


~  -  j  ( 1  +  erf  5) 


In  this  equation  the  central  velocity  line  where  —  =  —  ,  is  defined 

—  2 

U1 

on  the  x-axis  where  y  =  0  implying  that  ?  =  0  .  A  suitable  choice 
of  two  velocities  will  make  it  possible  to  calculate  the  similarity 
parameter  <r  .  The  following  method  is  us  d  in  this  study  to  find  the 
first  approximate  value  of  cr  . 

h 

Let  Y  denote  the  \  coordinate  of  the  points  where  —  =  a  . 

a  hj 

h  is  the  reading  of  the  pressure  manometer  which  corresponds  to 

the  dynamic  head,  and  subscript  1  refers  to  the  free  stream. 

h 

Similarly  corresponds  to  —  =  3  . 


Equation  — -  =  j  '1  +  erf  £) 

U1 


gives 


erf  ?  -  2  -  -  1  ---  2 


Ai  any  station  x  ,  ~~  a  gives  a  value  of  £  and  ~~ 

hi  °  hi 


gives  a  value  of  . 
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By  definition 


As  <r  is  assumed  constant  across  the  mixing  region,  the  similarity 
parameter  might  be  obtained  by 


e 

a 


AC 


Y 

a 

cr  — 


AY 

'"x 


Therefore 


Cft  ar.d  Cg  are  derived  tueoreticaily  from  the  equation  while  Y^ 
and  V,,  are  read  directlv  from  the  measurements.  This  method 

b 

makes  it  possible  to  determine  <r  from  the  empirical  data.  In  this 
thesis  analysis  is  based  on  cr  =  0.  9  ,  3  '  0.  t  .  Correspondingly, 
CQ  9  *  1-  *5  and  C^  t  =  -0.345  .  Yq  ^  and  YQ  ^  were  read  at 
each  station.  However,  this  approximate  value  of  a  did  not  yield 
the  best  agreement  between  theory  and  measurement.  It  became 
necessary  to  make  few  adjustments  to  improve  the  agreements.  The 
velocity  profile  based  on  the  measurements  were  compared  with  the 
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theoretical  profile  given  by  Eq.  (3-15)  in  Fig.  23  .  Results  of  the 
comparison  were  in  agreement,  which  implied  that  the  velocity 
profile  in  the  separation  region  ha;  the  same  characteristic  as  the 
half-jet  mixing  layer.  However,  it  was  discovered  that  the  simi¬ 
larity  parameter  <r  increases  with  distance  rather  than  remaining 
constant  for  all  value  of  x  .  <r  was  plotted  against  x  on  a  logarith¬ 
mic  paper  as  shown  in  Pig.  24  .  This  plot  showed  that  most  of  the 
points  were  scattered  around  a  straight  line,  except  those  very  near 
the  hill  and  those  approaching  the  reattachment  point.  A  slope  of 
approximate  0.  5  was  found  from  the  plot  which  meant  that  <r 
increases  witnvx  .  This  contradicts  the  assumption  leading  to 
Equation  (3-15)  in  which  cr  was  assumed  to  be  independent  of  x  . 

Despite  this  discrepancy,  if  cr  is  considered  as  a  pure 
empirical  factor  without  regard  to  its  physical  meaning,  Eq.  (3-15) 
still  gives  a  successful  description  of  the  velocity  distribution  along 
the  separation  region  behind  a  two-dimensional  hill.  However, 
further  insight  into  the  character  of  the  similarity  parameter  cr 
can  be  seen  by  considering  the  effects  of  the  initial  boundary  layer 
thickness  and  curvature.  With  the  consideration  of  the  initial  bound- 
ar\  layer  thickness  as  implied  in  Eq.  (3-17)  ,  Kirk's  approximation 
of  a  displaced  asymptotic  profile  will  modify  cr  into  cr 

x 

O . 

it  ait*  — ) 

x 
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The  displacement  of  the  virtual  origin  x  can  be  obtained  from 
Fig.  17  .  This  can  be  done  by  extending  the  straight  line  that 
represents  the  linear  spread  of  Zone  I  width  upstream.  The  inter  ¬ 
section  of  this  line  with  the  x-axis  is  then  the  virtual  origin  where 
it  is  assumed  the  half-jet  starts  to  develop.  The  distance  from  the 
crest  of  the  hill  model  to  the  virtual  origin  was  measured  and  was 
found  to  be  xq  =  20  inches.  Fig.  25  shows  the  similarity  para¬ 

meter  <r  obtained  by  this  method.  Although  a  is  still  not  a 
constant,  its  variation  with  distance  is  much  smaller  and  a  constant 
value  is  approached. 

The  variation  of  cr  with  x  might  be  attributed  to  the 

o  b 

curvature  of  the  separated  flow  which  is  due  to  the  transverse 
pressure  gradient.  This  effect  is  implied  in  Eq  (3-20)  . 

cr* 

"  \  1  -  mQ 

.  n  u  1 
where  Q  =  -  .  -r=- 

H 

9y 

j(l*  erf  ?) 


•  .  <lu 

it  Yields  — — 

<)  V 


_V_  t>-p 


'  n  x 


from  — 
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then 


Q  = 


u  xvr 

u(  R  a 


(5-6) 


It  is  obvious  that  the  value  of  Q  varies  with  y  .  It  was  decided  that 
the  mean  value  of  Q  could  represent  the  characteristic  of  the 
whole  region  across  the  separated  flow.  This  value  was  obtained  as 
follows: 


Q 

mean 
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V  i  Q  dy 
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where  y  and  ?  denote  the  value  of  y  and  ?  respectively 

V  ,  V  V/  , 

where  —•  =  0,9  ,  so  is  y_  .  and  .  The  integration  was 

h  0,1  0.  l  b 

carried  out  graphically  from  t  u.  345  to  ?q  9  "  1  -15 

it  yields 


C 

mean 


(5-7) 


Kor  /.er<»  curvature  the  similarity  parameter  <r 12  was  used 
The  value  of  m  can  be  then  obtained  from  Eq.  (  i-21)  .  Sawyer  used 
m  ;  5  in  his  analysis  ot  a  t w o-dumrisional  jet.  However,  the  value 
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of  m  found  i*~  tnis  study  is  not  a  constant  as  it  varies  from  5  to  6 
as  shown  in  Fig.  26  .  This  indicates  that  a  successful  description 
of  the  effect  of  curvature  on  the  turbulent  mixing  mechanism  is 
necessary. 


H.  Velocity  Distribution  in  Separated  Flow 

The  primary  purpose  of  this  thesis  is  to  investigate  the 
velocity  distribution  in  the  separated  flow  behind  a  two-dimensional 
hill.  The  development  of  mean  velocity  distributions  behind  the  hill 
can  be  seen  in  Fig  12  .  The  velocity  near  the  floor  is  not  plotted 
out  because  it  is  difficult  to  measure  the  reverse  velocity  in  the 
bubble  with  a  pitot  tube.  The  theoretical  velocity  profile  for  a  half¬ 
jet  is  given  by  Eq.  (3-15)  ,  which  is  used  to  compare  with  the  mea¬ 
sured  velocity  profiles  in  the  separated  flow.  From  measurements, 
free  stream  velocity  and  the  x-axis  can  be  determined  as  described 
previously  in  this  chapter.  A  suitable  value  for  the  similarity 
parameter  <j  is  then  chosen  and  the  measured  velocity  profiles  are 

plotted  out  non-dimensionallv  in  Fig.  ?3  with  7  against  ?  . 

U1 

Theoretical  velocity  profile  is  also  plotted  in  the  same  figure. 
Deviations  are  observed  in  the  lower  part  ot  Zone  II  .  This  deviation 
becomes  larger  near  the  floor  which  may  be  due  to  the  presence  of 
■;  .'-op:*'  nibble  and  the  teiiurt  oi  pitot  tube  to  measure  velocity 
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correctly  in  the  bubble.  However,  satisfactory  agreemen,  is  four.cl 
betweer  the  half-jet  theory  and  the  measurements  over  most,  of 
separation  region,  especially  in  Zone  I  . 
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Chapter  VI 

SUMMARY  AND  SUGGESTIONS 

A.  Summary 

Velocity  profiles  in  a  turbulent  separated  flow  of  an  incom¬ 
pressible  fluid  were  investigated  using  the  concept  of  free  turbulence. 
The  flow  separated  at  the  crest  of  a  wedge  shaped  model  hill  and 
reattached  itself  to  the  floor  downstream  of  the  hill  forming  a 
separation  bubble.  Gortler's  half-jet  theory  to  the  first  order 
approximation  is  employed  in  the  analysis  of  the  separated  flow. 
Analysis  is  simplified  by  neglecting  the  initial  boundary  layer  thick¬ 
ness  and  by  assuming  that  there  is  no  pressure  gradient  in  both  the 
main  stream  and  transverse  direction.  Satisfactory  agreements 
between  theory  and  experiments  were  louna.  The  following  results 
were  obtained: 

t.  The  velocity  profiles  in  the  separated  flow  over  most 
region  are  successfully  described  bv  the  half-jet  theory 
except  those  near  the  hill,  reattachment  point  and  the 
vortex  region  in  the  bubble. 


Initially  the  j»*t  boundary  and  jet  axis  are  curved  and 
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then  approach  to  straight  lines.  The  width  between  the 
two  lines  varies  linearly  with  the  distance. 

3.  A  similarity  parameter  cr  is  used  to  relate  the  measure 
ments  and  theoretical  values.  The  parameter  is  not  a 
constant  but  is  proportional  to  the  half- power  of  the 
distance  from  the  crest  of  the  model  hill. 

4.  The  curvature  of  the  separated  flow  is  induced  by  the 
action  of  the  pressure  gradient  across  the  mixing  region 
in  the  transverse  direction.  This  transverse  pressure 
gradient  is  large  only  in  a  short  region  immediately 
behind  the  hill.  Further  downstream  the  pressure 
gradient  becomes  so  small  that  the  curvature  of  the 
separated  flow  can  be  neglected. 

Finally,  the  similarity  parameter  is  modified  by  Kirk's 
approximation  of  displaced  origin  to  include  the  effect  of  the  initial 
boundary  layer  thickness.  The  effect  of  curvature  on  the  similarity 
parameter  is  approximated  by  the  first  order  theory  of  Sawyer. 
However,  the  results  indicate  that  a  better  interpretation  of  the 
effect  of  curvature  and  a  better  definition  of  the  radius  of  curvature 


are  needed. 
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B.  Suggestions  for  Further  Study 

It  is  recommended  that  further  studies  be  made: 

1.  to  measure  the  velocity  inside  the  separation  bubble. 
Since  the  pitot  tube  fails  to  give  a  reliable  measurement 
of  the  reverse  velocity  in  the  bubble,  an  improved 
instrument  will  be  needed  for  this  purpose. 

2.  to  find  the  streamlines  in  the  flow  field  behind  the  hill. 
These  were  not  obtained  in  this  study  because  no  mea¬ 
surement  was  made  to  read  the  reverse  velocity  in  the 
bubble. 

3.  to  determine  the  radius  of  curvature  of  the  separated 
How  in  a  more  logical  way.  In  this  study,  the  radius  of 
curvature  is  assumed  constant  across  the  separated  flow 
and  this  is  a  very  rough  assumption. 

4.  to  investigate  the  entrainment  of  the  fluid  from  the 
bubble  by  the  separated  flow  and  to  find  out  its  effect 
on  the  velocity  profile  in  the  separation  region. 

5.  to  include  in  the  analysis  of  the  separated  flow  the 
pressure  gradients  ir.  both  the  x-direetion  and  the 
v-direel  ion 
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6.  to  make  a  systematic  study  of  the  effects  of  the  shape 
and  of  the  size  of  the  hill  on  the  development  of  the 
velocity  profiles  in  the  separated  flow. 

7.  to  study  the  effect  of  the  standing  eddies  appearing  in 
front  of  the  hill  on  the  separated  flow. 
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FiG.  9  CALIBRATION  CURVE  OF  HOT  WIRE 


Equipment  used  in  measuring  reattachment  point 
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FIG  15  INTRINSIC  COORDINATES  AND  REFERENCE  COORDINATES, 


SPREAO  OF  THICKNESS  OF  ZONE 
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Fia  18  VARIATION  OF  STATIC  PRESSURE  ALONG  CONSTANT  VELOCITY 

RATIO  LINE  FOR  Urn  «  45  (p* 


VARIATION  OF  STATIC  PRESSURE  AT  CONSTANT  HEIGHTS  ABOVE  THE  FLOOR  WITH  UM*45fp* 
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DEVELOPMENT  OF  STATIC  PRESSURE  BEHIND  THE  HILL  WITH  U.«45fpv 


turbulent  components  behind  the  wedge  hill 


(in.) 


FIG  23-1  COMPARISON  OF  THEORY  AND  VELOCITY  MEASUREMENTS  FOR  U*  *  45  fpt 
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FIG  23-2  COMPARISON  OF  THEORY  AND  VELOCITY  MEASUREMENTS 


FIG  23-3  COMPARISON  OF  THEORY  AND  VELOCITY  MEASUREMENTS  FOR  U-,  *  45  fps 


FIG  24  VARIATION  OF  SIMILARITY  PARAMETER  a  WITH  DISTANCE 
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FIG.  26  VARIATION  OF  m  WITH 


